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SUMMARY ^ 

A numerical analysis of the eff ect of a circular fillet on the thermal 
radiation exchange between planar surfaces of zero conductivity intersecting at a 
sharp angle is presented. Temperature variation along the surface for a constant 
heat transfer and local heat-transfer variation for a constant wall temperature 
is obtained. Surface emissivity, which is assumed to be independent of wavelength 
and temperature, and fillet size are varied in the analysis. The "hot spot" for 
the unfilleted structure is found at the intersection of the planar surfaces but 
is moved to the point of juncture in the case of a filleted structure and is 
reduced in value. For planar surfaces that intersect at a right angle, the temper- 
ature and local heat-transfer distributions for various fillet sizes are compared 
with the corresponding distributions of the unfilleted surface. Exact results are 
compared also with the temperature and heat-transfer distributions obtained from 
an approximation to the radiation flux function which is based on the use of a 
shape factor. For a wide range of values for the emissivity and fillet size, the 
agreement between approximate and exact results is within 5 percent, n < ,r 

ffumorL 


INTRODUCTION 


Interest in radiative cooling of space vehicles, solar energy collectors, 
and thermal radiation shields has brought about a need to increase the fund of 
engineering knowledge through predictions of the characteristics of basic geometric 
configurations and through application of available mathematical techniques to the 
analysis of the fundamental equations. Heaslet and Lomax contributed to that fund 
by examining the radiative heat transfer between conducting planar surfaces that 
intersect at a sharp angle (ref. l) . These same authors also studied the radia- 
tive heat transfer for infinite, nonconducting shells with circular-arc sections 
(ref. 2). The purpose of this paper is to extend the results of these two recently 
published analyses by wedding the two configurations employed in their work. Thus, 
planar surfaces are joined by a concave surface of arbitrary radius, that is, 
joined by a fillet. The presence of the fillet at the juncture of the planar sur- 
faces affects the magnitude and distribution of radiation flux, .exact solutions 
for which are obtained from a machine code that incorporates the Li ouvi lie -Neumann 


method of solving a Fredholm integral equation of the second kind. Simple, closed- 
form solutions for certain limiting conditions provide a partial check on numerical 
results . 

The emphasis in this analysis will he placed on, first, the influence of the 
fillet on the radiative flux distribution corresponding to the unfilleted structure 
and, second, the importance of surface emissivity. Effects of thermal conduction 
will not he treated here, thus eliminating nonlinear terms in the governing inte- 
gral equation. The configuration, two-dimensional surfaces intersecting at right 
angles, suffices to illustrate the alteration of the flux distribution due to 
changes in the primary parameter, fillet size. As in references 1 and 2, the sur- 
face properties of absorptivity and reflectivity are assumed to be independent of 
surface temperature, and wavelength of radiation. The emissivity of the surface 
is assumed equal to the absorptivity, less than unity and constant at all temper- 
atures and at all wavelengths. By definition, then, the surface is gray. In 
addition, it is assumed that emission and reflection are diffuse and only the 
concave side of the configuration radiates. 

Two physical problems are formulated and solved for a filleted structure. 

In the first example the surface is held at a constant temperature and the varia- 
tion in local heat transfer along the surface is found. In the second example a 
constant heat transfer is provided and the variation in temperature along the sur- 
face is found. The solutions to these problems are compared with corresponding 
solutions for an unfilleted structure (some of which do not appear in ref. l) to 
demonstrate the influence of the fillet. 

To compare solutions corresponding to different fillet sizes, the total 
length of the combined surfaces is held constant. Fillet size, therefore, can be 
given simply as the ratio of fillet length to total length. When this ratio is 
1, the present results for local heat transfer and temperature variation are 
identical with results published in reference 2, which treats radiative transfer 
for infinite shells with circular-arc sections. For a length ratio of 0 the local 
heat-transfer solution without conduction is the same as that given in reference 1, 
but the second example, temperature variation with constant heat transfer for 
adjoining plates, is given here. 


NOTATION 


a see figure 1 

B(s) total radiation flux function (energy per unit area, time) 
c distance between field and source points 

f(s) shape factor defined in equation ( 31 ) 

G(s) unknown function in equation (12) 

H(s) flux of incident radiation defined in equation (3) 
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kernel function defined in equation (l) 

dimensionless length of straight-line segment on surface 
length of straight-line segment on surface cross section 
total length of surface cross section, 2S + 2L 
local heat-transfer function defined in equation (7) 


dimensionless radius of fillet. 


S + L 


S + L 


*R 


ratio of total fillet arc length to total surface length, 
radius of circular arc 

surface coordinate of point divided by S + L 
half-length of circular-arc segment on surface cross section, 

*+ 

temperature function 

Cartesian coordinate divided by S + L 
Cartesian coordinate divided by S + L 
absorptivity of gray surface 

dy 

slope angle of surface at field point , arc tan — — 

dx 

difference in left and right value of derivative of a function 

emissivity of gray surface 

angular coordinate of point on surface 

scale constant determined by equation (l8) when evaluated at s = 0 
parameter defined by equation ( 13 ) or (l4) 
reflectivity of gray surface 

angle measured from normal to surface at field point to line joining 

field and source points, arc sin (counterclockwise, +; 
clockwise, -) 
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Subscripts 


0 constant value of function 

R maximum value of function used as reference value 

1 approximate value of a function 


Superscripts 


source of radiation; unprimed coordinate refers to location of field point, 
that is, coordinate of receiver 


FORMULATION OF PROBLEM 


Consider a concave cylindrical surface of infinite width, a cross section of 
which is shown in figure 1, that is symmetric about a generatrix of the surface 
midway between the ends. Length along a cross section of the surface is measured 
from the midline and in a direction perpendicular to the generators of the surface. 
The characteristic dimension for this problem is the half-length, Lp /2, of the 
surface cross section and it is used to make all quantities involving length 
dimensionless. Since the surface is two-dimensional, an infinitesimal element of 
area of unit width is represented by ds or ds T where s' is the dimensionless 
surface coordinate that locates the source line element and s locates the field 
line element at which radiation is incident. 

The thermal radiation field (energy/unit time) incident at a field point, s, 
lying on the concave side of the surface is calculated by summing the total radia- 
tion emitted at all source points, s T , also lying on the concave side of the sur- 
face. Diffusely radiating elements are assumed to satisfy Lambert ! s cosine law. 
This means that when the normals to the elements, ds and ds', at the field point 
and source point, respectively, are projected onto the line joining the points, 
the flux received by each element is directly proportional to the area of the 
projected elements. In the idealized problem to be treated here, the normals, 

N and N l , are coplanar, affording a simple expression for the angle factor (see 
ref . 3 ) * 


— I sin cp( s, s' ) 

2 I I 

The incremental angle factor, or kernel function times the source element, ds r , 
determines the fraction of energy emitted and reflected at ds T that falls 
directly on ds and is given as 


k 


1 

2 


d( sin cp) 


= K( s, s T )ds 


a) 



where cp is the angle measured (counterclockwise, +: clockwise, -) from the 
normal to the surface at the field point to the line joining the source and field 
points. Equation (l) defines the kernel function, K(s, s T ) • The mathematical 
form and properties of the kernel that are pertinent to the present problem, apart 
from the matter of machine coding, will be set forth in the next section. 

Each element of the surface emits radiation according to the Stefan-Boltzmann 
law adjusted by a factor, e, the emissivity of the surface. An energy-flux balance 
for an element ds* on the concave side of the surface is written as a sum of two 
terms: the emitted radiation, eoT 4 (s ! ), and the reflected radiation, pH(s T ), where 

H(s ! ) is the energy-flux incident at s T and p is the reflectivity of the 
surface. The total flux, B(s ! ), ensuing from the source point is 


B(s*) = eaT 4 ( s’) + pH(s f ) 


( 2 ) 


where the function H(s T ) is 


H(s’) 


B(s)K(s, s’)ds 


L t 


(3) 


an integral of the total flux arriving at s ! summed over the entire length of 
the surface. The source and field points are interchangeable, so that the integral 
equation for the total flux leaving ds is 


B(s) 


ecrT 4 (s) + p / b(s*)K(s, s')ds' 

J L T 


( 4 ) 


For the gray surface considered here the coefficients of reflectivity, p, 
absorptivity, a, and emissivity, e, satisfy the relations 

1 = P + a (5) 

e = a (6) 

The latter equation is Kirchhoff*s law under the condition of temperature 
equilibrium. The heat-transfer function is defined as 

Q(s) = B(s) - H(s) ( 7 ) 

which, combined with equations (2), (5), and- (6) > gives 

B(s) = cT 4 (s) — Q(s) (8) 
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Inserting equation (8) into ( 4 ) casts the inhomogeneous integral equation ( 4 ) 
into a form convenient to formulate the two problems mentioned earlier. The basic 
equation which serves as a starting point for the numerical work is 


€ctT 4 (s) - Q(s) = J [eaT 4 (s’) - pQ(s')]k(s, s')ds' 


and the two problems are stated next in dimensionless form. 

If a constant temperature, T 0 , is maintained everywhere on the surface, then 
equation ( 9 ) becomes 


e + (l - e) P [ ^ - s 4 " — “ — 1 K(s, s T )ds 


lyp [ 0 T 0 1 - e 


If, instead, a constant heat transfer, Qq, is provided everywhere , then equation 
(9) takes the form 


^ 4 (s) 1 , p crT 4 (s T ) _ 1 - € 


'LT l ^ 


+ / gSLILg < - — - — I K(s, s T )d.s r 


In the formal mathematical solution of equation (10) or (ll) or both, one 
considers the Fredholm integral equation 


G(s) = 1 + A / G( s t )k( s , s T )bs 

J T.m 


where, for the constant temperature case, 


€ £ 2 CJT 0 


A = 1 - e 


and in the case of constant heat transfer 


gT 4 (s) 1 - £ 
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A = 1 


It can be seen from equations (12) and (l4) that the gradient of the temperature 
field is independent of the emissivity, whereas the gradient of the heat-transfer 
field is altered by the presence of the factor (l - e)/e 2 in equation (13)* 


Equations (10) and (ll), or equivalently, equation (12) with (13) or (l4), 
are solved numerically for the configuration shown in figure 2- It consists of 
planes of length L whose extensions intersect at right angles which are joined 
by a concave surface of length 2S, a section of which is a circular arc. The 
parameter, r g , is defined as 


S 

r s " S + L 


(15) 


and will be varied in the analysis along with the emissivity. In order to compare 
different configurations, the total length of the surface, Lp, where 


L T = 2S + 2L 


(16) 


is held constant. Since dimensions involving length have been made dimensionless 
by dividing by the surface half-length, Lp/2, the limits of integration in the 
surface coordinate system become +1 and -1. 

Before examination of the composite kernel function in the next section, it 
is noted that equations (10) and (ll) have a simple analytic form in the special 
case that 


r s = 1 


that is, when the plane surface length, L, is zero. The configuration then is an 
infinite shell with a circular-arc section for which the solutions to the two 
problems stated above are given in reference 2. They are 


= /cos S I + ^ tan | cos cos f s') (17) 

j-j-rn ^ \ O 0 O/ \ O / 


for the constant shell temperature and 


aT*(s) _ 1 , 7 t 

Q 0 = e + 8 





(18) 


for the uniform heat transfer. The symbol R, the radius of the circular arc, is 
related to the fillet parameter, r s , by the expression 


R = 


2r s L p 

Jt 


(19) 
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From this identity the angular coordinate, 0, employed in reference 2 and 
equations (l7) and (l8) here, is related to the surface coordinate and fillet 
parameter by 


0 = f- 


( 20 ) 


for j s | < r s • 

Exact solutions to equations (10) and (ll) are obtained from equations (12), 
(13) j and (l4) with the aid of a Fortran program for the IBM 7090. 


KERNEL FUNCTION 


The kernel function that appears in the integral equation governing the 
radiation flux function is defined as the absolute value of the gradient of the 
angle factor. For any two-dimensional, continuous configuration with piecewise 
continuous curvature, the angle factor is 

|sin cp(x, y; x' , y')l = i fe P ^ (y - y’) s in 4 (2l) 

J(x - x* ) 2 + (y - y* ) 2 


where 3 = arc tan dy/dx. A transf ormation from Cartesian coordinates to surface 
coordinates (see fig. 2) is convenient when dealing with a composite surface. 

For three choices of location of field and source points, equation (2l) assumes a 
compact form. For example, when both points lie on the circular arc, the angle 
factor is 


1 cos — — |s - s' I 

2 8r s ' 

and is employed in reference 2. When both points lie together on either straight- 
line segment, the angle factor is unity; whereas if the field point is on one 
straight part and the source point is on the other which is the case considered 
in reference 1, the angle factor (eq. (2l)) can be written as 
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The different forms of the kernel function corresponding to the three angle 
factors enumerated above as veil as the angle factors which result from the other 
possible locations of field and source points are computed using equation (l) and 
are listed in figure 3* A perspective drawing, figure A, depicts the approximate 
behavior of K(s, s 1 ) for r s = 0.3 and exposes a property of the composite kernel 
function that is an aid to checking the final numerical results. It is evident 
that the gradient of K(s, s f ) is discountinuous at several places. At the 
juncture of the concave and plane surface, the radius of curvature jumps discon- 
tinuously from a finite value to infinity. Hence, the derivative of the kernel 
is discontinuous at the juncture also, namely, at s* = r s and s = r . (it is 
also discontinuous at s = s T .) 

The discontinuity in the kernel gradient causes a break in the slope of the 
curves that represent Q(s)/crT 0 4 and crT 4 (s)/Q 0 at s = r g . At the juncture point 
the derivative of either function takes on two values which are determined accord- 
ingly as s = r s is approached from the left or from the right. The slope differ- 
ence formula for each function can be found, first, by writing down the total 
radiation flux, B(s), as seen from points on either side and very close to the 
juncture, then by differentiating with respect to the field coordinate, s, and 
finally, by passing to the limit, r g . Thus, let si and s 2 be points to the left 
and to the right of r s and very close to it. The difference in slope of the 
functions, Q(s)/T q 4 and T 4 (s)/Q q , can be obtained most simply from the dimension- 
less, canonical form of the basic integral equation, (12), with equations ( 13 ) and 
(lA), in the place of the dimensional flux function B(s) mentioned above. 
Expressions for G(s) to the left and right of r g are 


G( s ) — 1 + A 


r p- r s p Si 

J g(s’)Kjy*(si, s*)ds* + J g(s’)Kj(si ,s’)<3.s' 

. “1 " r R 


+ / G(s t )Kj(sij s^d-S* + j G(s T )Kjy(si, s f )ds 


'si 


( 22 ) 


and 




G(s 2 ) = 1+Al / G( s t )Kjjj*(s 2 , s')d.s' + / G(s , )%j(s 2 , s T )ds f 




J G(s*)KiIl( s 2 ^ s')ds' 


(23) 


The slope difference formula for G(r q ) is 


c)G(s) 


_ lim 5 g(si) lim 5g( s 2 ) 
s=r s s i"* r s Ssi s^+r s ^ S2 


(24) 
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and when equation (24) is combined with (22) and ( 23 ), the final form is written 
compactly as 


6 


3G(s) 


3s 


ls=r s 



G(s’) 


dKjv* 
^ 3s 


^mi) ds- 

3s / 


+ 





ds* 



^ K IV 

3s 



(25) 


Expressions for K(s, s T ) are given in figure 3* 

As can be seen from figure 4 the composite kernel function itself is 
continuous, is never negative, and is bounded everywhere except in the limit as 
r g goes to zero and s = s T = 0. Note that a factor of 1/2 is put into the pres- 
ent definition of the kernel which is not present in the functions displayed in 
references 1 and 2. The singularity that occurs when the arc length or fillet 
radius is zero and s = s T = 0 is illustrated in reference 1. The function, 
Kiii*(s, s * ) , for r s = 0, is not bounded throughout the entire region of integra- 
tion but, rather, increases indefinitely as one approaches the origin along an 
arbitrary line, s/s* m constant. However, the value of the function, G(s), at 
the inner corner of the adjoint plates is finite and is determined with the aid 
of a Li ouvi lie -Neumann expansion in powers of the parameter, A. The sum of this 
series is a simple formula which was found and exploited to great advantage in 
reference 1. For the two examples elaborated here G(0) yields 


Q(0) = e 
ctT q 4 1 + e 

and 

oT 4 (0) = 1 + e 
Q o “ e 


(26) 

(27) 


These formulas are of considerable value in tying down the numerical solution in 
a region where the computed data begin to scatter, that is, near s = s T =0 for 


Another special situation, not involving a singularity of the kernel, arises 
when r s = 1. Here K(s, s 1 ) reduces to Kj(s, s 1 )* the inversion of equation (12) 
is possible, as shown in reference 2, and formulas (17) and (l8) result. 
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PRESENTATION OF RESULTS 


Two problems of physical interest have been formulated so that they emerge 
directly from equation (4) when a transformation, equation (8), is made. The 
transformation casts (4) into a form containing the two functions, Q(s) and a!T*(s), 
explicitly. When either function is regarded as a given constant of the problem, 
the basic equation for the other function reduces to equation (10) or (ll). It is 
only for convenience in presentation, however, that one function is regarded as a 
constant. The machine program developed to handle equation ( 9 ) and applied to (10) 
and (ll) is equally valid when an arbitrary boundary condition is given. To facil- 
itate the presentation of the numerical results that illustrate the effect of a 
fillet, however, it suffices to use a constant, either T 0 or Qq, as the prescribed 
function on the boundary. Solutions to the integral equations, (10) and (ll), 
discussed here may be found in figures 5 aa d 6, respectively. The emissivity and 
fillet size appear as parameters. 

Figure 5* a g^aph of Q(s)/crT 0 4 versus s, consists of three groups of 
curves: an upper group, middle group, and lower group, illustrating the level of 

variation for three representative values of the emissivity, namely, e = 0-9, 0. 5 , 
and 0 . 1 . Each curve within a group depicts the actual behavior of the function. 
Several values of fillet size are plotted, including the two limiting cases, 
r s = 0 and r s = 1.0. 

It is evident from figure 5 that, for a given surface temperature, T Q , the 
heat-transfer distribution, Q(s)/cxT 0 4 , is affected in two ways. First, for any 
value of emissivity the distribution for the unfilleted structure tends toward 
uniformity when a fillet is added. The spread between maximum and minimum heat- 
transfer values of any curve diminishes with increasing fillet size. The great- 
est spread occurs for the unfilleted structure whereas the circular arc (r g = l) 
sustains the least variation of heat transfer along its surface. Second, for a 
given fillet size the level of the variation of the heat transfer diminishes with 
decreasing emissivity. As one might expect, therefore, a fillet has the greatest 
effect for surfaces with high emissivities . The influence on surfaces with low 
emissivities is small. 

Figure 6 presents the temperature distribution on a filleted structure with 
constant heat transfer. The scale constant, ® Q , that appears on the figure is 
the value of the function oT 4 (o)/Q o at s = 0 for r s = 1 and for a given e 
(eq. (l8)); it is presented as figure 6(a). As € decreases 0 O increases and 
the curve is arbitrarily terminated at e = 0.01 as a practical limit. The tem- 
perature function aT 4 (s)/Q 0 - 0 O is given in figure 6(b) along with three 
dashed curves. The latter curves are typical results based upon the use of an 
approximation which introduces a shape factor. A discussion of the approximation 
is given in the next section. 

It can be seen from the results In figure 6(b) that for a given heat 
transfer, Q , to the surface, the temperature distribution, cfT 4 (s)/Q 0 , is affected 
only by the fillet size and not by the surface emissivity, a consequence of the 
transformation that was found, namely, equation (l4) . The maximum temperature 
occurs at the corner of the unfilleted structure but this temperature is reduced 
and occurs at the juncture of the two surfaces when a fillet is introduced. For 
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a fixed value of emissivity the spread "between maximum and minimum temperatures 
diminishes with increasing fillet size. The greatest spread occurs for the 
unfilleted structure whereas the circular arc, as before, sustains the least var- 
iation in temperature along its surface. In this case, however, the gradient 
along the surface is independent of emissivity. It should be noted though, that 
the temperature variation along the surface becomes less and less significant as 
the emissivity is reduced since the spread in temperatures becomes a smaller and 
smaller fraction of the maximum temperature of the surface. Thus, as with the 
previous example, the influence of the fillet is reduced with decreasing 
emissivity. 

For the range of parameters explored in this presentation only small values 
of e cause convergence difficulty. A method of accelerating convergence was 
developed in reference 1 and is contained in the machine program used here. Two 
other parameters that have not been mentioned so far but appear in the machine 
program are the step size and the desired accuracy of the computed result. A 
total of 51 equi spaced field points on the surface at which a value of a function 
is calculated are used. This is a step size of 0.02 since the total surface 
length, Lt, is set equal to 1 in the machine code but equals 2 in the formulation 
of the problem and the presentation of computed results. A finer step size, 0.01, 
did not alter any numerical results, to four significant figures, except when 
very near the origin for small fillet sizes. Some variation in the third signif- 
icant figure then is observed but little difficulty is experienced in fairing 
curves because the value at the origin is known. The iteration procedure, for 
all computed results, was continued until the difference in successive iterates 
was less than 0.0001. 

The value of either function, Q(s)/crT 0 4 or ctT 4 (s)/Q 0 , at the juncture of the 
fillet and plane surfaces is obtained by extrapolation of the curves on either 
side of the juncture point. This is necessary inasmuch as the juncture point is 
not one of the equispaced field points. Except for fillet sizes less than triple 
the step size, that is, 0.06, this is satisfactory because there are then three 
points through which the curve can be constructed. However, it is not necessary 
to use a finer step to obtain the curve corresponding to fillets less than 0.0 6 
because the value of the function is known at the point of symmetry, s = 0, and 
the value at the junction can be obtained by extrapolating only the solution from 
the right. The ordinates of the intersection of this extrapolated curve and the 
ordinate, s = r s , is the desired value. These two values provided an upper and 
lower bound on the cosine function which corresponds to the fillet portion of the 
surface. The slope of this function is zero at s = 0 and there is enough infor- 
mation, in most instances, to estimate adequately the curve between s = 0 and 
s = r s . In the two figures, 5 and 6(b), curves for r s = 0.01 and 0.05 are drawn 
in this way. 

The code that was developed by H. Lomax for equation (9) with kernel 
function, equation (2l), differs from that discussed in reference 1 in the manner 
in which the conf iguration and composite kernel are introduced. Also it was 
unnecessary to incorporate into the program the analytic formula to replace the 
corner singularity. However, the iterative scheme and convergence tests are the 
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same and the reader may consult this reference concerning these two aspects of 
the numerical method as well as a discussion of the technique by which the rate 
of convergence can be accelerated. 

It is interesting to consider an alternate presentation of the solutions to 
equations (10) and (ll) which define two distinct physical problems. Instead, 
one might begin with a mathematical problem defined by equation (12), the 
Fredholm integral equation of the second kind with symmetric kernel. The kernel 
is characterized by the parameter, r g , as before. The solutions to this equation 
can be represented by a two-parameter family of curves, G(s), the parameters 
being A and r s instead of e and r s . An example of this family is given by 
figure 6(b) when e = 1.0. This group of curves represents, in the mathematical 
sense, the variation of G(s) with s for A = 1.0 and several values of the kernel 
parameter, r s . The envelope of the solutions shown in this figure is defined by 
the curves corresponding to r s = 0 and r s = 1.0. As the parameter A is 
reduced from 1 to 0, the size of the envelope as well as the level of the entire 
group of curves is also reduced. In the limit that A = 0, the envelope coa- 
lesces to the single curve, G(s) = 1.0. The transformations, equations (13) and 
(l4), lead from the mathematical problem as stated above back to the two physical 
problems which have been discussed previously and whose solutions are given in 
figures 5 and 6. 


APPROXIMATION USING SHAPE FACTOR 


An approximate solution, Gi(s), of the canonical form of the radiation 
integral equation, equation (12), is obtained under the assumption that the vari- 
ation of G(s') with s' can be ignored with respect to the variation of K(s, s' 
with s' . It is possible, therefore, to replace G(s’) by G(s) and the resultant 
expression is 


Gi( s) 

where f(s), called the shape factor, is 


1 - Af( s ) 


(28) 


or alternately 


f(s) 



s* )ds’ 


(29) 


and 


f(s) = 1 - ± 
2 


sin cp(s, -l) - sin cp(s, +l) 


s > 0 


(30) 



(31) 


See, for example, figure 1. Equation (30) is especially convenient to graphical 
determination of the shape factor . 
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For comparison with the exact numerical solution of equations (10) and (ll), 
the corresponding approximate solutions are 


Qi(s) _ [el- f(s)] 
cT 0 4 1 - (1 - e)f( s ) 

and 

oTi 4 (s) _ 1 f(s) 

Q 0 el- f(s) 


(32) 


(33) 


As a measure of accuracy of the approximation, the quantities 


Q( s) - Qi( s ) 

% 

and 


(34) 


a T 4 ( s ) - crT 1 4 (s) 
0T R 4 


(35) 


ve re computed for a representative number of points and for the same values of 
the two parameters illustrated in figures 5 and 6 (b). For these values of the 
two parameters, the absolute value of (34) is found to be less than 0.09 and the 
absolute value of ( 35 ) is found to be less than 0.05 everywhere. 

In both physical problems considered, the greatest disagreement between 
exact and approximate solutions occurs for a fillet size of about 0.1. The 
dependence of expressions (34) and (35) on emissivity, in contrast to their 
dependence on fillet size, differs for the two cases. For the constant tempera- 
ture case, expression (34) increases with decreasing emissivity, being 0.006 for 
e = 0.9, whereas it is 0.09 for e = 0.1. For the constant heat-transfer case, 
expression (35) decreases with decreasing emissivity. An emissivity, e = 0.9, 
yields the value 0.05 for expression ( 35 ) but when € = 0 . 1 , the same expression 
yields a value of 0.009- These results suggest also that in the mathematical 
sense the approximation, Gi(s), most accurately represents the exact value G(s), 
given by equation ( 12 ), for low values of A where 0 < A < 1 . Approximate 
results, shown as dashed lines on figure 6 (b), are presented only for the con- 
stant heat-transfer case because this case represents a significant departure of 
the approximate solution from the exact solution and at the same time the largest 
variation in the function crT 4 (s)/Q 0 . The constant temperature example for high 
values of e in which the variation of Q(s)/T 0 4 is significant is well 
represented by the formula, equation (34). 

At the corner, s = 0, for all values of emissivity and for r s = 0, the 

approximation G-i(s) is seen to yield exact results; that is, equations ( 30 ) and 

( 32 ), and ( 30 ) and ( 33 ) lead to equations ( 26 ) and ( 27 ), respectively, since 

f(0) = 1/2 for r g = 0. This explains, in part, why a maximum is observed in 

expressions (34) and ( 35 ) when they are regarded as continuous functions of r . 

s 


i4 



Without the corner formulas, (2 6) and (27), the available numerical data used in 
equation (29) and, hence, the approximate solutions, (32) and (33), do not suggest 
such a maximum. The approximation appears to break down whereas it is actually 
exact at the corner. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 24, 1962. 
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Figure 3.- Kernel function in surface coordinates listed according to relative 
location of field (s) and source (s') points. 







